contact impurity potential [22]) 

•s.=xM-£)E 



where the spectral function is defined as 
sgnco ImE n 



A(k,oo) 



71 (co-^ k -ReEj2 + (ImEj 



2 ' 



(5) 



(6) 



and /(co) the Fermi function. At low T, 3//3co converges 
to a delta function, and only the self-energy around E p is 
important. We assume that the 3D coupling in the cuprates 
is sufficient, so that localization effects can be neglected in 
the low density limit under consideration. 
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Fig. 4. (a) Frequency dependence of the normal state self- 
energy Im£ n for various values of l/u. (b) Impurity induced 
resistivity at T = 100K, as a function of 1 /u for 8 = 0.18,0.25. 

Fig. 4(a) shows the CD dependence of Im£ n for 8 = 
0.25, n. = 0.005, and various values of u. The surprising ob- 
servation is that even though the overall shape of the curves 
is different, in the vicinity of E F , they are very similar. In 
particular for l/u = -0.1,0.7 the values at t p are essen- 
tially the same. The vHs creates some type of an inversion 
point for the frequency dependent scattering rate. This 
similarity reflects itself in the impurity induced dc resis- 
tivity which we plot as a function of 1 /u, for n. = 0.01, T « 
100K, and 8 = 0.25,0.18 in Fig. 4(b). The values be- 
tween l/u = -0.5,0.7 vary by less than 15%, even though 
the T suppression within this range changes by a factor 
of two [see Fig. 3(c)]. The magnitude also agrees well 
with reported values for Zn and Ni (« 20 - 30pI2cm/at%) 
[6]. This shows that the proximity to the vHs, causing the 
strong co dependence of E n is capable to severely violate the 
conventional proportionality of residual resistivity and T 
suppression found for a constant DOS. The large resistiv- 
ity caused by Ni together with its mild effect on T as com- 
pared to Zn, might be explainable by this effect. 

In conclusion, we have shown that a simple BCS model 
with a d 2 _ 2 OP, supplemented by the experimental quasi- 

x y 

particle dispersion, and non-magnetic impurity scattering 
provides a consistent description of the experimentally ob- 
served effects of Zn and Ni. The presence of the vHs 



slightly below ^ F is necessary to resolve the puzzle given 
by the violation of the proportionality between residual re- 
sistivity and T suppression for the two types of impurities. 
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as l/u = —0.1, i.e., \u\ « 1.5eV ss W, the bandwidth, is suf- 
ficient to explain the values reported for Zn [4]. This as- 
signment immediately leads to the prediction that the Zn 
induced A^ res should decrease significantly for underdoped 
samples, e.g., for 8 = 0.13, N res (n ( - = 0.002) « 0. In Fig. 
1(b), we plot the n. dependence of N for selected val- 

l II 

ues of 1 /u, and 8 = 0.25. The expected relation A^ res ~ n ( ' 
holds for l/u = —0.1, whereas for l/u> 0.6, iV is neg- 

' ' ' ' res ° 

ligibly small up to concentrations of n. = 0.05. The latter 
behavior was reported for Ni doped Y-0 ? , hence, we as- 
sign a value of l/u « 0.7 to this case. Note, that if sub- 
stitution occurs primarily on the planar Cu site, n. « 1.5x 
in YBa 2 (Cu t _T X ) 3 7 _ y (T = Zn.Ni). However, at larger 
x, substitution is also expected on the Cu(l) chain sites. 
Fig. 1(c) confirms the point made earlier, that the occur- 
rence of resonant scattering is now linked to the condition 
|c(C0 = 0)| <C 1. For 8 = 0.25, |c(fi) = 0)| =0, i.e., resonance, 
occurs for 1 fu « -0.0425. 
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Fig. 2. The to dependence of the DOS for small n. = 0.001 near 
resonance (a,b), and for larger n ( . far from resonance (c,d). 

Figs. 2(a-d) show the frequency dependence of the 
DOS near e F . For very small n. [Figs. 2(a,b)], we obtain 
the characteristic hump at e provided the resonance con- 
dition |c((0 = 0)| <C 1 is satisfied. On the other hand, for 
weaker scattering l/u = 0.7 [Figs. 2(c,d)], N((t>) vanishes 
linearly at £ p for n j < 0.02, and for larger n., a sublinear 
dependence is found with a small N , which grows as 8 is 
reduced. From this, we predict that a small N ks could ap- 
pear in underdopedNi substituted cuprates for n. < 0.05. 

To calculate T , we solve the linearized gap equation in 
the presence of impurities 



y (cosk —cos A; ) 2 

^5 



(4) 



The pair potential is chosen as V kk , = Vr| k r| k , to generate 
a d x 2_ y2 OP and (3 c = (k B T c )~ l . Due to the non-trivial t, k , 
the normal state (A = 0) self-energy entering (4) becomes 
frequency dependent. Hence, T cannot be expressed in 



the standard AG form, and the Matsubara sum has to be 
evaluated numerically with a cut-off, typically a> n < 50 (« 
7.5eV). For the pure case, this results in an accuracy bet- 
ter than 10~ 4 compared to the value without cut-off. The 
coupling constant was chosen as V = 2.2, which corre- 
sponds to a near neighbour attraction of V/S = 0.275 ss 
41meV (see [18]). Without impurities, this leads to T = 
(111K,92K,65K) for 8 = (0.25,0.18,0.13). 
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Fig. 3. (a,b) The T suppression for potentials near and far 
from resonance, for 8 = 0.18,0.25. (c) l/u dependence of the 
T suppression rate for 5 — 0.13,0.18,0.25. (d) Frequency de- 
pendence of Im£ () on the imaginary axis at T for various 1/u. 



Figs. 3(a,b) show the n. dependence of T at 8 = 
0.25,0.18 for the values of l/u assigned to Zn and Ni. The 
shape of the curves is essentially identical to the AG form 
(see for example [20]). The T suppression rate is about 
two times as fast for l/u = -0.1 (Zn) than for l/u = 0.7 
(Ni). Furthermore, it strongly increases for lower 8, in 
qualitative agreement with experiments on YBa 2 Cu 3 ? 
[21]. The critical concentrations at l/u = —0.1 (n c « 0.06 
for 8 = 0.25, n L i « 0.04 for 8 = 0.18) also agree well with 
experimental values when taking into account the planar 
site preference of Zn, as discussed above. Fig. 3(c) illus- 
trates the substantial increase in dt c /dn. [t = T c /T c0 ,T cQ = 
T c (n i = 0)] upon lowering 8. In analogy to Af [Fig. 1(a)], 
the peak shifts from negative to positive l/u as 8 is low- 
ered. Fig. 3(d) shows the frequency dependence of the self- 
energy E Q on the imaginary axis for various values of 1 /u 
[Re{E (r )} = 0, and Im{E 3 (7;)} = 0]. The large variation 
of dt c /dn. between 1 fu = — 0. 1 , 0.7 originates in the big dif- 
ferences of E Q particularly at intermediate co n . 

Finally, we need to check whether the u values assigned 
to Zn and Ni can produce a quasiparticle damping ImE , 
which could explain the similarity in the residual resistiv- 
ity caused by the two. The normal state dc conductivity is 
given by the Kubo formula (vertex corrections vanish for a 



3 



formed more easily than Cu +m . (iii) The vHs also results 
in a strong frequency dependence of the normal state quasi- 
particle lifetime. This might explain why T suppression 
rate and increase of residual resistivity are not simply re- 
lated as observed for Zn andNi. (iv) Both, N , andofr jdx 
are strongly doping dependent, the latter becoming much 
larger in the underdoped regime away from the vHs. (v) 
Overall, good quantitative agreement with experiments is 
obtained with only one adjustable parameter, u. 

We model the impurities by a short-range potential 
V(r) = u8(r - r.), and apply the self-consistent t-matrix 
approximation [8, 16]. This approach has recently been 
shown to yield accurate results in the dilute impurity limit 
[17]. Working in particle-hole space, the t-matrix T satis- 
fies the following Lippmann-Schwinger equation (quanti- 
ties with a hat represent matrices) 



e vHs * s published to date) respectively, e vHs the position of 



vHs 



T((0) = u 



a 3 + S 3 f(a))-£^(k,(o) 



(1) 



Here we introduced the propagator 

o5a + A k S 1+ ! k a 3 



f(k,co) 



(2) 



where co : 



: & ~ E 0> A k = A k + E 1 ' = + E 3' the C l Uasi " 

particle energy. The self-energy (like all other matrices) 
is expanded in terms of the identity and Pauli matrices 
a Q • • • a 3 as E = EyS,-. and is given by E = n.T, n i the impu- 
rity concentration, and A the volume. For the d x2 _ y2 OP, the 
off-diagonal self-energy vanishes by symmetry, and from 
(1), the non-zero components are 



n ; G n ( .(l/w — G 3 ) 

E ° = (l/«-G 3 )2-G2 ; L i= { l/ u -G 3 f-Gi 



(3) 



where G ( .(co) = A -1 £ k g.(k,co). Previous approaches [5, 
8] assumed particle-hole symmetry in the normal DOS, 
which has the consequence that E 3 = is a self-consistent 
solution of (3). Here, we relax this condition, and take into 
account the full structure of the t-matrix. Unfortunately, 
this also means that all k-sums have to be done numerically 
when iterating Eqs. (3) for a self-consistent solution. A 
d x d k-space grid with d = 4000 for calculations of spec- 
tral quantities, and d = 200 to calculate T proved to be suf- 
ficient to obtain good convergence. 

For the quasiparticle dispersion we use a 

tight-binding fit to ARPES data on Bi-2212 with 
real space hopping matrix elements [r ,---,r 5 ] = 
[0.879,-1,0.28,-0.087,0.094,0.087], (t Q on-site, t x nn, t, 
nnn hopping, ... ) [18, 19]. All energies are measured in 
units of |f j | = 0. 149eV. The value of t Q corresponds to a hole 
doping of 8 = 0. 17. In the calculations, we consider doping 
levels of 8 = 0.25,0.18,0.13 (t Q = 1.0,0.9,0.8) which are 
representative for Y-0 ? (e yHs ss -16meV [15]), Bi-2212 
(e vHs « -31meV [14]), and Y-O g6 (no precise value of 



the vHs with respect to e A rigid band picture is assumed 
to be an acceptable approximation in this small doping in- 
terval. In the calculations of spectral quantities, we set 
A = 0.2 (30 meV). The only free parameter in the theory 
is then the potential strength u, since the impurity concen- 
tration n. is also fixed by experiment. We shall search for 
values of u to consistently explain the experiments for Zn 
and Ni. 

The effect of the non-trivial DOS is rather dramatic, in 
particular if a vHs is close to £ p , as in the present case: (i) 
The s-wave scattering phase shift 8 Q acquires a strong fre- 
quency dependence already in the normal state, which also 
reflects itself in the scattering rate, hence resistivity. Fur- 
thermore, the dependence on u is highly non-trivial, (ii) In 
the self-energy Eqs. (3), the cotangentc = cot8 (being/re- 
quency independent in case of a constant normal DOS) that 
is usually used to parameterize the scattering strength, is 
replaced by the (now frequency dependent) quantity c(co) = 
- 1 ju + G 3 (co) . This leads to a strong sensitivity of the su- 
perconducting DOS and T on both u and the chemical po- 
tential. Resonant scattering is usually observed for \c\ <C 1 
(corresponding to \u\ « °° for constant DOS). In the present 
case, this translates to the condition c(co = 0) 1 , allowing 
for resonant scattering even if \u\ <C °°. We shall illustrate 
these points in the following. 
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Fig. 1. (a) The residual DOS N , as a function of 1/u and dif- 
ferent doping levels, (b) N res as a function of impurity concen- 
tration n j for various values of 1 ju. (c) The real and imaginary 
part of c(co) for values of l/u near resonance. 

Fig. 1(a) shows A^ res = N(co = 0) = -7t" 1 Im{G (co) + 
G 3 (co)}| (B=0 , as a function of the potential strength u for 
fixed n. = 0.002. A broad peak is observed, the center of 
which shifts from negative (1/u « —0.05 for 8 = 0.25) to 
positive (1/u « 0.1 for 8 = 0.13) values as the hole doping 
is lowered. The broadness of the peaks shows that the sig- 
nature of resonant scattering, i.e., alargeA r res for very small 
n., occurs in a rather extended interval on the 1 /u axis. Us- 
ing 8 = 0.25 to model Y-0 ? , we find that a value as small 
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We consider the effect of non-magnetic impurities in a d x2 _ y2 superconductor with e p close to a van Hove 
singularity. It is shown that the non-trivial density of states (DOS) allows for resonant scattering already 
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The effect of non-magnetic impurities on the properties 
of a superconductor (SC) can provide useful information 
about the symmetry of its order parameter (OP). Ander- 
son's theorem [1] states that non-magnetic impurities nei- 
ther affect the transition temperature T , nor the density of 
states (DOS) N((0) of a BCS SC with an isotropic OP. How- 
ever, if the OP is anisotropic, possibly exhibiting nodes 
(such as e.g. a d x2 _ y2 OP), non-magnetic impurities may 
lead to a large residual DOS N vss = N((£> = 0) at the Fermi 
energy e F , and even to a complete suppression of T [2]. 

In the high-T copper-oxides, Zn impurities are partic- 
ularly harmful. For YBa 2 (Cu 1 j Zn x ) 3 1 for instance, 
where Zn substitutes primarily on planar Cu sites [3], a 
rapid suppression of T (T = at x ss 0.08 — 0.1), as well as 

a jV. es growing like x 1 / 2 is observed [4]. On the other hand, 
nominally magnetic Ni +2 impurities (also believed to pri- 
marily substitute for planar Cu) have a much milder effect: 
while suppressing T at a rate 1/2 - 1/3 of Zn, they do not 
lead to a significant iV res up to x « 0.05 [4]. 

An appealing and simple interpretation of these re- 
sults is possible in terms of a d x2y2 OP affected by non- 
magnetic scattering, assuming that Zn acts as a strong, in 
fact resonant scatterer [5] . However, one problem with this 
model comes from the fact that the growth in the normal 
state resistivity dpjdx just above T is quite similar for Ni 
and Zn [6]. Given that, their different behavior in the SC 
state is puzzling, since in standard models (for magnetic 
and non-magnetic impurities) [7,8], the T suppression rate 
dT c /dx, and the residual resistivity are both determined by 
the same parameter, the normal state scattering rate. Also, 
it is not clear why Zn should scatter resonantly, since this 
is obtained only for scattering strength \u\ -» °°. 

Reports about localized magnetic moments on Cu sites 
next to a Zn [9] further complicate the issue. The size of 
the moment depends on hole doping, and seems to be large 
in underdoped cuprates (« 0.9 — 1.4jj, b ) [9, 10], whereas 
in optimally doped YBa 2 Cu 3 7 y (Y-0 ? ) it is small [9] or 
maybe absent [11]. Therefore, at least in Y-0 7 , it seems 



justified to model Zn as a pure potential scatterer. In under- 
doped samples, additional magnetic pair breaking might 
occur. Note however, that even then it is not obvious that 
Abrikosov-Gorkov (AG) theory [7] applies, due to (i) the 
existence of short-range antiferromagnetic (AF) correla- 
tions, and (ii), the fact that the 'impurity' spins are proba- 
bly not totally immobile, but part of the same SC spin fluid 
(Cu spin 1/2). For a strongly correlated system like the t-J 
model for instance, it was shown [ 1 2] , that a localized spin 
1/2 impurity only has a small influence on magnetic and 
pairing properties when its coupling to the mobile spins /. 
is the same or similar to the coupling J among the latter. 
This should be satisfied here. There is also the possibility 
that Zn might act on the pair potential itself as proposed re- 
cently [13] in the context of an AF spin fluctuation model. 
However, the presence of localized Cu moments nearest 
neighbour to Zn suggests that AF correlations are proba- 
bly unaffected or maybe even enhanced in the vicinity of 
the impurity site, since this region corresponds to a locally 
underdoped phase, i.e., is closer to AF. 

As one possible route to resolve the discrepancies be- 
tween theory and experiment described above, we explore 
the influence of a realistic quasiparticle dispersion on the 
effect of non-magnetic impurity scattering in a model BCS 
superconductor with an assumed d x2 _ 2 OP chosen as A k = 
Ar| k , with r| k = (cosk x — cosfe y )/2. Our main findings 
are: (i) the existence of a van Hove singularity (vHs) just 
below £p (suggested by angular-resolved photoemission 
[ARPES] [14, 15]) leads to a strong violation of particle- 
hole symmetry, and allows for the occurrence of resonant 
scattering, at realistic potential strengths \u\ « 1 — 2eV. 
(ii) The sign of u is crucial. For hole concentrations cor- 
responding to Y-0 ? , we find that resonant scattering oc- 
curs only for negative u (in electron notation). A large 
negative u is expected for Zn impurities, since an inert d- 
shell strongly repels holes, i.e., attracts electrons. This is 
also consistent with the observed localized Cu moments. 
Ni, however, should create a weaker attractive potential 
for holes, i.e., positive u, since a Ni +m oxidation state is 
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